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SUMMARY 

NASA-Lewis and the Navy .iointly developed, based on a previous Navy code, 
an engine cycle program capable of simulating the design and off-design per- 
formance of arbitrary turbine engines. Boeing, under a NASA-Lewis contract, 
developed a computer code which, when used in conjunction with the cycle code, 
could predict the weight of the engines. Thus, the Navy/NASA Engine Program 
(NNEP) along with WATE-2 (Weight Analysis of Turbine Engines) determines the 
dimensions and weight of each major component in the engine. The output from 
these codes originally included a very crude and indecipherable computer 
drawing of the flowpath representation of the designed turbine engine. This 
proved to be of little use. It was desired to add a graphics subroutine to 
the code to enable the engineer to visualize the designed engine with more 
clarity by producing an overall view of the designed engine for output on a 
graphics device using IBM-370 graphics subroutines. In addition, with the 
engine drawn on a graphics screen, the program would allow for the interactive 
user to make changes to the inputs to WATE-2 for the engine to be redrawn and 
reweighed. These improvements would allow better use of the WATE-2 code in 
conjunction with NNEP. 



INTRODUCTION 


NNEP (Ref. 1) is a very general cycle analysis code that can assemble an 
arbitrary matrix of fans, turbines, ducts, shafts, etc. into a complete gas 
turbine engine and compute on- and off -design thermodynamic performance. 

WATE-2 (Ref. 2) is a preliminary design procedure for calculating engine 
weight using the component characteristics determined by NNEP. A main objec- 
tive of the WATE-2 program using NNEP is to get a "feel" for what the total 
designed engine looks like in terms of size, weight, and its flowpath. Previ- 
ously, the only output that could be examined as being the "total engine" was 
the flowpath configuration, shown in Figure 1, which was very hard to deci- 
pher. Development of a subprogram that incorporates the IBM-370 graphics 
package subroutines has now allowed the capability for a clearer output on a 
graphics device. The graphical output of the engine includes numbers on the 
components, a display of the total engine's weight, and an option of individ- 
ual data output of the individual components which includes type, number of 
stages, and weight of the component, as shown in Figure 2. 

DISCUSSION OF METHODS 

Interactive Mode 

In the past, if an engine was to be modified, the program running was 
terminated, input dataset was modified, and the program was rerun. Now avail- 
able is the direct interactive creation and modification of the engine's 
inputs without exiting from the program. The interactive mode allows modifi- 
cation of any or all inputs to WATE-2 for revisions to the engine. The inter- 
active mode is used in conjunction with the graphics subroutine (EGPLOT). 

Thus, the flag for graphics must be turned on in the weight code input 
(PL0T=T). The program will prompt the user on Unit 20 as to whether modifica- 
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tions are desired, thus Unit 20 must be identified to the computer as being 
the interactive terminal. Request of modifications will key the program to 
prompt the user for the input modifications. The new NAMELIST inputs entered 
on Unit 20 will over-write the previous inputs and proceed to perform the 
weight-code estimate. 

Additions to the subroutine WTEST in the code WATE-2 are shown in Appen- 
dix C (marked with an arrow at each line in the subprogram that is new or has 
been modified). 

Example Case 

Appendix D contains a photocopy from an actual terminal session for inter- 
actively modifying a component for a more reasonable-designed engine. Shown 
in Figure 3, component Number 6 (low pressure compressor) is off-set to the 
fan (No. 2). The beginning of Appendix D shows the original WATE inputs, with 
the inputs to be changed marked. The H/T ratio of the compressor is increased 
(DESVAL (3,6)) and the mean diameter of duct 4 increased (DESVAL (3,4)). 

Thus, the newly designed engine is shown in Figure 4 with the compressor 
lined-up with the fan. 

Non-Interactive Mode 

The non-interactive mode allows the user to perform the WATE-2 program but 
not make any changes to the inputs. The code is written so that it will run 
as in the past, i.e., without the graphics subroutine, unless the user inputs 
a flag for graphics. Graphics is performed by putting PL0T=T in the weight 
code input (&W). PL0T=F will turn off the flag and the graphics will not be 
used. When executing the program in the non-conversational mode (BACK) the 
program and the procedure definition have been developed to perform with or 
without graphics. If plots are desired, then PLOT must be set to "True," as 
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above, but also an input to the procedure definition, WATE2IT, must be 
BACK=YES. The output of the plots will be on the Zetal2 plotter. If BACK is 
set to "YES" but the program is not executed in background, the program will 
still output plots on the Zetal2 plotter and will not prompt for any input. 
Description of Subprogram 

An explanation of the subprogram and subroutines is given below. A more 
thorough explanation of the graphics subroutines is given in Appendix A. 


GRAPHICS SUBROUTINES 


Line 102 

GRINIT 

Initializes graphics package 

Line 103 

CORNER 

Removes corner marks on PLOT 

Line 131 

SETFRM 

Sets no. of frames 

Line 148 

XAXIS 

Produces the x-axis 

Line 153 

YAXIS 

Produces the y-axis 

Lines 164, 217, 

GPLOT 

plots the given x and y array values 
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Line 186 

GARC 

Draws an arc or circle 

Lines 260, 275- 

NIMBER 

Converts a number to an array of 

277, 287 


Printable graphics characters lEBCDIC) 

Lines 261, 270, 

CHARS 

Positions printable characters 

278-282, 


(EBCDIC) anywhere on the graphics 

286, 288 


Display 

Line 263 

TITLE 

Puts a title on the PLOT 

Line 289 

DISPLA 

Displays the accumulated plot orders 


on the graphics device. 
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Program Description 
Lines 105-135 


This section of the program finds the maximum radius of the engine 
and scales the y-axis according to the maximum radius. The user is prompted, 
if in the conversation mode, for single or multiple frames and whether compo- 
nent output is desired on the plot. The x-axis is scaled to the engine's 
length and the x-axis is set according to the desired frame option. 

Lines 136-153 

Using the previous scaling information, the x- and y-axis are setup 
for the graphics output. 

Lines 155-170 

The centerline of the engine (radius=0) is created for the plot. 

Lines 171-1^11 

A curve for an arbitrary inlet is created using a generated arc. The 
slope of the line for the inner radius of the second flowpath component is 
matched and the arc generated from that point to the centerline. This part of 
the program is not necessarily needed but could help in visualizing the engine. 

Lines 192-258 

This section draws the basic components of the engine to be output on 
the graphics device. "Ducts" are drawn with a cosine function to create a 
"smoother" looking engine to allow for misalignment of components. This is 
the main section of the program. The components are drawn in the order of the 
flowpath (JFLOW) using output summary data from WATE-2. 

Lines 259-271 

The component numbers are put on the components to help in under- 
standing the structure of the engine drawn. Also, the title of the engine is 
put on the plot. 
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Lines P72-290 


If component output on the plot was requested, this section of the 
program will be performed. It converts data into a printable graphics charac- 
ter array that is put on the plots. 

Lines 291-301 

The total engine weight is put on the plot and then the final graph- 
ics output is displayed on the designated output graphics device. 


CONCLUSIONS 

The modifications to WATE-2 described in this report can help visualize 
needed changes to the engine's inputs and reduce time spent modifying the 
inputs. The user can interactively change inputs such as H/T ratio, aspect 
ratios of the blades, number of compressor or turbine stages, etc. to redraw 
the engine and see the effect on engine dimensions, flowpath and weight. The 
user can thus be assured that the final engine design is reasonable. 

The logical structure developed for this code is adaptable to any other 
code in which the dimensions of the components are known. 

NNEP and WATE-2 computer codes are available by written request to L. H. 
Fishbach, c/o NASA Lewis Research Center, 21000 Brookpark Road, Mail Stop 
501-10, Cleveland, Ohio 44135. Requesters will be expected to input the 
necessary graphics subroutines for their system to replace the IBM 370 graph- 
ics subroutines indicated in this subprogram (EGPLOT). 
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Figure 2 "New" graphics displayed flowpath representation (full length 
scale; reduced). 
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Figure 3 Sample engine before interactive changes (NOTE: Component 6). 
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Figure 4 Sample engine after interactive modifications to inputs. 
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APPENDIX A 

Subroutine EGPLOT User-Guide for Interactive Mode 

The following user-guide pertains only to the IBM-370 computer system at 
NASA Lewis Research Center. However, this guide will show the steps taken 
when using the interactive mode for all systems. Use of the graphics devices 
mentioned and the set-up of the system's graphics package may differ from 
system to system. 

The computer program is written so that it will run as in the past, i.e., 
without the graphics subroutine, unless you input a flag for graphics. This 
is done by putting PL0T=T in the &W input for the weight-code. PL0T=F will 
turn off the flag and the graphics will not be used. 

When PLOT is set "TRUE" the graphics subroutine, EGPLOT, will be 
referenced. If you do not want all of the plots displayed, you will still 
have the option of graphics or no graphics. You will be prompted as to 
whether you sti''! want the engine plotted, i.e., "DO YOU WANT A GRAPHICS 
PICTURE? YES = Y; NO = N." If you do not, then it will go on to the next 
scale factor (or terminate) without performing the graphics routine. If you 
request an engine graphics display, then it will use the graphics routine. 

The graphics picture of the engine can be arranged in a few ways. You 
will be prompted by the program for the desired arrangement. The first option 
is whether you want multiple frames or not, i.e., "DO YOU WISH MULTIPLE 
FRAMES? YES = Y; NO = N." If you do not want multiple frames then you will 
get just one frame which has the components out of proportion with the x-axis 
condensed. The multiple frames gives a scaled picture with the engine in 
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proportion. Usually the single frame picture of the engine is sufficient for 
analysis of the engine, but multiple frames will give a better looking engine, 
especially on a Zeta or Calcomp plot. The user will also have the option as 
to whether component output is to be printed on the picture or not, i.e., "DO 
YOU WANT COMPONENT OUTPUT? YES = Y; NO = N." If the component output is 
desired, the component numbers along with its weight, type, and number of 
stages will be put on the engine plot. The total weight of the engine is 
still put on the plot even if the component output is not. 

If a graphics device for output has not been previously defined, then it 
must be defined before the display of the engine can be performed. The system 
will prompt you for a device. Enter the desired output device. This device 
will stay the same for all plots until you "abend" or do a GSUnCH and 
redefine the output device. 

Upon completion of the displayed engine, you will be prompted as to 
whether changes to the weight-code's inputs are desired, i.e., "DO YOU WISH TO 
MAKE CHANGES TO THE INPUT? YES = Y; NO = N." If using a graphics terminal for 
execution of the program and also for the graph cs display, t^is question may 
come before initializing the graphics mode. If nothing is entered before 
displaying the picture, just enter the appropriate response after exiting from 
the graphics mode. If the bell has been turned off the terminal, backspace 
one and then enter the response.. The prog>^am will go on ’-'ith the 
weight-estimation program, if no changes are requested. If modifications to 
the inputs are requested, then the program will prompt you for the new 
NAMELIST inputs, i.e., "ENTER DESIRED CHANGES FROM TERMINAL, I.E., &W DESVAL 
fl,6') = *****, &END." Any number of changes can be made to the NAMELIST 
inputs. The new NAMELIST inputs will over-write the previous inputs and again 
carry-out the weight-estimation program. 
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When using offline graphics printing, such as the Calcomp or the Zeta 
plotters, GTERM must be entered (interactive mode) to initialize the batch 
task for offline printing to obtain the desired plots. 
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APPENDIX B 


PROGRAM LISTING FOR SUBPROGRAM EGPLOT 


SUBROUTINE EGPLOT 


PURPOSE 


TO MAKE GRAPHICS PLOT OF THE TOTAL ENGINE SHOWING INDIVIDUAL 
COMPONENTS ALONG THE FLOWPATH CJFLOW). 


DESCRIPTION 


THE OVERALL VIEW OF THE DESIGNED ENGINE IS PRODUCED FOR OUTPUT ON 
A GRAPHICS DEVICE USING I3M-370 GRAPHICS SUBROUTINES. IN THE 
INTERACTIVE MODE ( BACK= ’ NO ' ► PLOT=T ) THE USER IS PROMPTED FOR SIZE 
OF THE PLOT ( PROPORTIONAL LENGTH OR SINGLE FRAME) AND A REQUEST 
FOR COMPONENT DATA OUTPUT ON THE PLOT. IN RUNNING THE PROGRAM IN- 
BACKGROUNDCBACK) ,BACK=»YE5» MUST BE AN INPUT TO THE PROCDEF AND 
PLOT=T SET. THE COMPONENTS ARE DRAW AS IT FOLLOWS THE FLOW PATH. 
USAGE 

CALL EGPLOTCENGLEN^WATENG) 

CALLING ROUTINES 


WTEST THE WEIGHT ESTIMATION ROUTINE 

REQUIRED SUBROUTINES 


>>>>>IBM-370 GRAPHICS PACKAGE SUBROUTINES<<<<<< 

GARC - DRAWS AN ARC OR A CIRCLE 

GRINIT - INITIALIZES THE GRAPHICS PACKAGE 

CORNER - REMOVES CORNER MARKS FROM PLOT 

5ETFRM - SETS THE NUMBER OF FRAMES FOR THE PLOT 

XAXI5 - CREATES A NONSTANDARD X-AXIS 

YAXIS - CREATES A NONSTANDARD Y-AXIS 

GPLOT - PLOTS DATA FROM ARRAYS 

NUMBER - CONVERTS A NUMBER TO CHARACTER DATA FORM (EBCDIC) 

CHARS - PLACES CHARACTER DATA ANYWHERE ON THE PLOT 

TITLE - PRINTS AN AXIS OR PLOT TITLE 

DISPLA - DISPLAYS THE ACCUMULATED PLOT ORDERS 
ON THE GRAPHICS DEVICE SELECTED 


MODIFICATION HISTORY 


DATE ID ANALYST DESCRIPTION 

MO/DA/YR IDEN7 NAME DESCRIPTION OF CHANGES 

AUTHOR/LANGUAGE/DATE 

ROBERT CORBAN-COOP (U. OF CINCINNATI) /FORTRAN IV / SEPT 21,1981 
GLOSSARY 

NAME ORIGIN USAGE DESCRIPTION 

IWMEC /WMECH/ I CONTROL INFORMATION 

WATE /WMECH/ I WEIGHT OF EACH COMPONENT 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
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ALENG 

TLENG 

RI 

RO 

ENGLEN 

WATENG 

VARS 

IVARS 

ITICKS 

NFRAM 

AMAX 

SKIPIT 


/WMECH/ I 
/WMECH/ I 
/WMECH/ I 
/WMECH/ I 
ARG I 
ARG I 
L 
L 
L 
L 
L 
1 


ACTUAL LENGTH OF EACH COMPONENT 

ACCUMULATED LENGTH TO END OF COMPONENT 

RADIUS INNER INLET, OUTLET EACH STATION 

RADIUS OUTER INLET, OUTLET EACH STATION 

TOTAL ENGINE LENGTH 

TOTAL ENGINE WEIGHT 

ARRAY OF DATA FOR DESIRED AXIS 

ARRAY OF CONTROL INFO. FOR LINE DRAWING 

NO. OF TICK MARKS ON THE AXIS 

NO. OF FRAMES FOR GRAPHICS OUTPUT 

MAXIMUM RADIUS OF ALL COMPONENTS 

BACKGROUND FLAG ( BACK=YES ) , SKIP PROMPTS 


SUBROUTINE EGPLOTC ENGLEN , WATENG) 

REALxa DATINP,DAT0UT,WTF,T0PRE5,T0TEMP,FAR,C0RFL0,VMACH,STATP,- 
lERROR,TOL,TOLT,TOLTT,DEPV,DTOL,PERPF,RPMNT,TMTEMP,TMPRES,DATOUM,- 
2DATMAC,DATALT,DE5LIM,TNPRES,TNTEMP,CNRFL0,C0RFLM,WNTF,DAT0UN,FARN,- 
3DANINP, DEBUG, DEPQ,SEL AST, DD,TOLOPT,M 

K COMMON BLOCKS X 

COMMON /DBL/ DATINP(15,60),DATOUT(9,60),WTF(<fO),TOPRES(^0),TOTEMP- 
l(AO),FAR('fO),CORFLOC90),VMACH(AO),STATP(AO),ERROR(40),TOL,TOLT- 
2,TOLTT,DEPV(20) ,DTOL(20) ,PERPF(20) 

COMMON /SNGL/ JM1,JM2,JP1,JP2,JCX,LOCTBL(9,60),JCOMP(70),IUAY,NIT,- 
1ITABC70),JCONF(60,A),JTYPE(60),JFLOW(70),IDEDAP(15),KKINDS<14,25),- 
2NCOMP,NOSTAT,NITER,NFINIS,NPASS,JCC,NTBL,NCTS,JCIND(20),JCDEP(20),- 
3JCVIND(20) , JCVDEP(20),KDTYPC20) ,IDONE(60) 

COMMON /WMECH/ IWMEC(7,60),WATE(60),ALENG(60),TLENG(40),RK2,40),- 
lRO(2,40),DESVALa7,60),DSHAF(5),RPMT(60),IWT,IPLT,IERR,ISII,ISIO,- 
2I0UTCD,NSTAG(60) 

COMMON /SKIP/ SKIPIT 

DIMENSION VARS(IO) ,X(30),IVARS(10) ,Y(30),SY(5),IS(60),ANGLE(3) 
DIMENSION DATA(7),WEIGHT(4),XS(60),YSC60) 

DATA DATA/4HC0MP,4H# ,4HTYPE,4H N,4HSTAG,4HE ,4H WT / 

DATA WEIGHT/4HTOTA,4HL EN,4HGINE,4H WT=/ 

DATA INLT, YES,C0MP/4HINLT,2X4HY / 

DOUBLE PRECISION PL0T3,PL0T4 
LOGICAL SKIPIT 
IWMEC(1,1)=INLT, 

CALL GRINIT ' 

CALL CORNER(l) 

FINDS MAXIMUM RADIUS OF ENGINE CAMAX) 

AMAX=0.0 

MAIN=0 

PHI=3. 1415927 
DO 20 I=l»2 
DO 10 J=l,40 

IF(R0(I, J) .GT. AMAX) AMAX=RO(I,J) 

10 CONTINUE 
20 CONTINUE 


PROMPTS FOR TYPE OF PLOT 

IST0RE=AMAX/5. 

YMAX=ISTOREX5.+5.0 

N=YMAX/5.+2. 


61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 
^5 
86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
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NFRAM=1 

IF CSKIPIT) GO TO 30 
WRITE(20,500) 

READ(20,600) FRAME 
WRITE(20,700) 

READ(20,600) COMP 
IF(FRAME.HE.YES) GO TO 30 
ITICKS=ENGLEN/5.+6 
IFRAM=N/7.X9. 

NFRAN=ITICK5/IFRAM 
CALL SETFRM(NFRAM) 

GO TO AO 

30 ITICKS = EMGLEN/10 .+3 
I0DD=ITICK5/2K2 

IFdOOD.NE. ITICKS)ITICKS = ITICKS+1 

SET-UP OF X AND Y AXIS 

AO VARS(1)=7. 

VARS(2) = 9.)<NFRAM 
VARS(3)=0. 

VARS(A)=-10. 

VARS(5)=(ITICKS-1)X10. 

IF(FRAME.EQ.YES) VARS(5) = ( ITICKS-2)i<5 . 

XMAX=VAR5(5) 

VAR5(6)=ITICK5 

VARS(7)=2. 

CALL XAXIS(0,5r0.5,VARS) 

VAR5(2)=7. 

VAR5(3)=90. 

VARS(5)=YMAX 

VARS(6)=N 

CALL YAXI5(0.5,0.5,VARS) 

CREATION OF THE CENTERLINE FOR THE ENGINE ( ) 

N=ITICKS-2 

X(1)=0.0 

X(2)=8.0 

Y(1)=0.0 

Y(2)=0,0 

IVARS(1)=3 

IVARS(2)=2 

DO 50 1=1, N 

CALL GPLOTCX,Y,IVARS) 

X(l)=XC2)+5. 

J = I/2 

IF(I.EQ.2?^J) X(2)=X(2) + 13. 

IF(I.NE.2HJ) X(2)=X(2)+9. 

50 CONTINUE 

CREATE AN ARBITRARY INLET USING A GENERATED ARC 

I=JFL0W(2) 

IUP1=JC0NF(I,1) 

IDN1=JC0NF(I,3) 

YO=RI(l, IUPi)K7.0/(YMAX+lfl. ) 

A = (RI(2, IDNl)-Ria,IUPl))X7.0/(YMAX+10. ) 

B = 9«NFRAM/(XMAX+10 . )?^ALENG(2) 

M=A/B 
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150 
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159 
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DEGREE=PHI/2.-ATAN(M) 181 

IF(DEGREE.GT. PHI/2.) DEGREE=PHI/2 . 182 

RAD=YO/SIN(DEGREE) 188 

CENTX=SQRT(RADXX2-YOX«2)+0.5+10.X9KNFRAM/(XMAX+10.) 18A 

CENTY=10.x7.0/(YMAX+10.)+0.5 185 

ANGLE(1)=0.0 186 

ANGLE(2)=90 . t ( PHI/2 . -DEGREE) 180 . /PHI ) 187 

ANGLE(3)=180.-ANGLE(2) 188 

CALL GARC(CENTX,CENTY,RAD,ANGLE,2) 189 

190 

COMPONENTS ARE DRAWN FROM SUMMARY DATA (RI AND RO) 191 

192 

DO 170 11=1,60 193 

I=JFLOW(II) 194 

NC=JTYPE(I) 195 

IF(I.EQ.O) GO TO 180 196 

IUPl=JCONF(I,l) 197 

IDN1=JC0NF(I,3) 198 

X(1)=TLENG(IDH1)-ALENG(1) • 199 

X(2)=X(1) 200 

SY(2)=R0(1,IUP1) 201 

SY(3)=R0(2,IDH1) 202 

IF(NC.EQ.8) SY(2)=SY(3) 203 

IF(MAIN.EQ.l .AND. NC.EQ.9) GO TO 80 204 

IF(NC.EQ.9) MAIN=1 205 

SY(1)=RI(1,IUP1) 206 

SY(4)=RI(2,IDN1) 207 

GO TO 90 208 

80 AREA=DAT0UT(5,I) 209 

SY(1)=RI(2,IST0R) 210 

SY(4)=SQRT(SY(3)XSY<3)-AREA/PHI) 211 

90 AINC=ALENG(I) 212 

IST0R=IDN1 213 

IVARS(2)=2 214 

N=4 215 

IF(NC.EQ.2) H=1 216 

DO 100 IP=l,N 217 

Y(1)=SY(IP) 218 

Y(2)=SY(IP+1) 219 

IF(IP.EQ.4) Y(2)=SY(1) 220 

CALL GPLOT(X,Y,IVARS) 221 

X(1)=X(2) 222 

X(2)=X(2)+AINC 223 

AINC=AINC-ALENG(I) 224 

100 CONTINUE 225 

IF(NC.NE.2) GO TO 160 226 

227 

CREATE A SMOOTH DUCT USING THE COSINE FUNCTION 228 

229 

JJ=JFLOW(II+l) 230 

110 CONTINUE 231 

1F(ALEN6(JJ).NE.0.0 .OR. JJ.EQ.60 .OR. JTYPE(JJ) .EQ.8) GO TO 120 232 

JJ=JJ+1 233 

GO TO 110 234 

120 IUP1=JC0NF(JJ,1) 235 

IUP2=JC0NF(JJ,2) 236 

SY(3)=RI(1,IUP1) 237 

IF(JTYPE(JJ).E0.9) SY(3)=RI(2,1DN1) 238 

SY(4)=R0(1,IUP1) 239 

1F(JTYPE(JJ).EQ.8) SY(3)=RI(1, IUP2) 240 
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IF(JTYPECJJ) .EQ.8) 5Y( ^ ) =R0 ( 1 , IUP2 ) 241 

DO 150 IJ=1,2 242 

DIFF=5Y( IJ+2)-SY(IJ> 243 

YO=SY( IJ)+DIFF/2. 244 

YMAX = ADS(DIFF/2. ) 245 

THEATA=0.C 246 

IF(DIFF.GT.O.O) THEATA=PHI 247 

DO 140 IK=1,30 248 

Y(IK)=YriAX>COS(THEATA)+YO 249 

IF( IK. EQ . 1 )GO TO 130 250 

X(IK)=X(IK-l)+ALENG(I)/29. 251 

130 THEATA=THEATA+0 . 1047198 252 

140 CONTINUE 253 

IVAR5(2)-30 254 

CALL GPLOT(X,Y, IVARS) 255 

IF(IJ.EQ.l) YSAVE=Y(15) 256 

150 CONTINUE 257 

C 258 

C LABEL COMPONENTS BY COMPONENT NUMBER (JFLOW) . 259 

C 260 

160 Yl=( 5Y( 1 )+5Y(2) )/2. 261 

IF(^K. EQ. 2) Yl=(Y(15)+YSAVE)/2. 262 

X1=TLENG( IDNl )-ALENG(I)/l .5 263 

IFCALENGCI) .LE.10.0 .AND. FRAME . NE . YES ) XI =X1-ALENG( I)X0 . 5 264 

IFCALENGf I) .EQ.0.0) 60 TO 170 265 

IF(X1 . GT . EHGLEN) XI =TL ENG( IDNl )~AL ENG( I )X0 . 95 266 

CALL NUMBER(1,I,2,0,NUMS) 267 

CALL CHARS(-2,NUMS,0.0,X1,Y1,10) 268 

170 CONTINUE 269 

180 CALL TITLE(1,60,15,IDEDAP) 270 

C 271 

C COMPONENT OUTPUT DATA IS PLACED ON THE PLOT, IF REQUESTED(C0MP=YE5) 272 

C - 273 

IF(COMP.NE.YES) GO TO 200 274 

X2=4.0 275 

YY=9.8 . 276 

CALL CHAR5(28,DATA, 0 . ,X2,9.8,15) 277 

DO 190 111=1,60 278 

I=JFLOW(III) 279 

IF(I.EQ.O) GO TO 200 280 

YY=YY-.2 281 

CALL NUMBER(4,WATE(I) ,6,0,PLOT3) 282 

CALL NUMBER(1,I,2,0,PLOT1) 283 

CALL NUMBERCl ,NSTAG(I),2,0,PLOT2) 284 

CALL CHARS(2,PLOT1,0. ,4.2,YY,ld) 285 

CALL CHARS(4,IWMEC(1,I),0:,5.3,YY,10) 286 

CALL CHARS(2,PLOT2,0. ,6.5,YY,10) 287 

CALL CHARS(6,PLOT3,0. ,7.5,YY,10) 288 

190 CONTINUE 289 

C 290 

C WEIGHT OF TOTAL ENGINE IS PLACED ON THE PLOT; PLOT IS DISPLAYED 291 

C 292 

200 CALL CHAR5(16, WEIGHT, 0. ,0.2,8.5,15) 293 

CALL NUMBER(4,WATENG,6,0,PLOT4) 294 

CALL CHARS(6,PLOT4,0. ,2.5,8.5,15) 295 

CALL DISPLA(l) 296 

500 FORMAT (’ DO YOU WISH MULTIPLE FRAMES? YES=Y ; NO=N') 297 

600 F0RMATCA4) 298 

700 FORMAT (' DO YOU WANT COMPONENT OUTPUT? YES=Y ; NO=N') 299 

RETURN 300 

END 301 
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APPENDIX C 


PROGRAM LISTING FOR REVISED SUBPROGRAM WTEST 


2 

SUBROUTINE WTEST 3 

4 

5 

PUJ?POSE 6 

7. 

TO CONTROL THE CALLING OF SUBROUTINES WHICH WILL ESTIMATE THE 8 

WEIGHT AND LENGTH OF INDIVIDUAL COMPONENTS 9 

10 

DESCRIPTION 11 

- 12 

THE OVERALL LENGTH OF THE ENGINE IS CALCULATED BY PROCESSING THE 13 

ILENG ARRAY. ALL COMPONENTS EXCEPT DUCTS AND SHAFTS, THEN DUCTS. 14 

THE REMAINING COMPONENTS EXCEPT DUCTS AND SHAFTS ARE PROCESSED. 15 

THE DUCTS ARE PROCESSED AND FINIALY THE SHAFTS. 16 

A BUILT-IN ASSUMPTION IN THE DUCT ROUTINE IS THAT NO DUCT IS • 17 

CONNECTED TO ANOTHER DUCT I.E. THE DUCT SIZE IS DETERMINED BY THE 18 

ADJOINING COMPONENTS. 19 

THEN THE MAXIMUM RADIUS IS FOUND. THEN DEPENDING ON THE PRINT 20 

FLAG -lOUTCD- THE REQUIRED PRINTING IS DONE 21 

IF THE PLOT CODE FLAG -IPLT- IS TRUE ROUTINE EPL*T IS CALLED 22 

USAGE 23 

24 

CALL WTEST 25 

26 

CALLING ROUTINES . 27 

28 

FLOCAL- 29 

ZTOPZ - 30 

31 

REQUIRED SUBROUTINES 32 

33 

COMP -COMPRESSOR WEIGHT/LENGTH 34 

TURB -TURBINE 35 

SHAFT -SHAFT 36 

DUCTW -DUCT 37 

COMBWT -PRIMARY BURNER WEIGHT/LENGTH 38 

WTNOZ -NOZZLE WEIGHT/LENGTH 39 

WMIXR -MIXER 40 

WSPLT -SPLITTER 41 

EPLT -PRINTER/PLOTTER 42 

43 

MODIFICATION HISTORY 44 

45 

46 

DATE ID ANALYST DESCRIPTION 47 

- - 48 

MO/DA/YR IDENT NAME DESCRIPTION OF CHANGES 49 

9/21/81 ROBERT CORBAN INCORPORATE INTERACTIVE MODE AND GRAPHICS 50 

51 

AUTHOR/LANGUAGE/DATE 52 

53 

NORMAN PREWITT-BOEING COMPUTER SERV. /FORTRAN IV / OCT 10,1976 54 

55 

GLOSSARY 56 

57 

NAME ORIGIN USAGE DESCRIPTION 58 

— — - 59 


IWMEC /WMECH/ I CONTROL INFORMATION 


60 
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WATE /WMECH/ 0 

ALENG /WnECH/ 0 

TLENG /WMECH/ 0 

RI /1‘jr^ECH/ 0 

RO /WMECH/ 0 

DESVAL /WMECM/ I 

D5HAF /WnECH/ 0 

RPMT /WMECH/ I 

IWT /UfiECH/ I 

IPLT /WMECH/ I 

PLOT L 

I ERR /i'-riECH/ 0 

IS 10 /mECti/ 1 

1 51 1 /WMECH/ I 

lOUTCD /WilECH/ I 

ILENG /WriECH/ I 

IDID 0 

<^KIPIT L 

SUBROUTINE WTE5T 

REAL 5(8 DATINP,DAT0UT,WTF,T0PRE5,T0TEnP,FAR,C0RFL0,VMACH.STATP,ERR0- 
lR,TOL,TOLT,TOLTT,DEPV,DTOL,PERPF,RPr,NT,TnTEnP,TMPRE5,DATOUM,DATNAC- 

2, DATALTrDESLIM,TNPRES,TNTEMP,CNRFL0,C0RFLM,WNTF,DAT0UN,FARN,DANINP- 

3, DEBUG,DEPQ,5ELA5T,DD,T0L0PT 

COMMON BLOCKS ^ 

XKXKXKHMXXKKXXKXX 

COMMON /DBL/ DATINP(15,60),DATOUT(9,60),WTF(A0),TOPRES(40),TOTEMP(- 
1AO),FAR(AO),CORFLO(AO) ,VMACH(AO),STATP(AO),ERROR(^0),TOL,TOLT,TOLT- 
2T, DEPVCZO ) ,DTOLC20 ) ,PERPF(20) 

COMMON /SNGL/ JM1>JM2,JP1,JP2,JCX,LOCTBL(9,60),JCOMP(70)|IWAY,NIT,- 
1ITAB(70), JC0NF(60, A) , JTYPE(60) , JFLOWC70) ,IDEDAP(15) >KKINDS(1^,25),- 
2NC0MP , NOSTAT,NITER,NFINIS,NPAS5/JCC,NTBL,NCT5,JCIND(20),JCDEP(20),- 
3JCVINDC20), JCVDEPC20).KDTYP(20),IDONE(60) 

COMMON /DEFAUL/ DEFAUL ( 1 5 , 2 0 ) , ISCAL E( 3 ) , SCAL E( 6 ) 

COMMON /WMECH/ IWMEC(7,60)>WATE(6i0),AL':NG(60),TLENG(Ad),RI<2,A0),R- 
10(2,AO) , DESVAL Cl 7,60) ,D5HAF( 5), RPMTC 60), IWT, IPLT, lERR, 1511,1510, 10- 
2UTCD,N5TAG(60) 

COMMON /CONVER/ C0NVERC15) 

COMMON /NEPCPT/ DEBUG,DEPQ,SELA5T,DD,T0L0PT,NDSET,NPART5,I0PTP,NPA- 
1S50,NV0PT,NJ0PT,N0PT 

COMMON /TERMON/ TNPRE5(AO),TNTEMPCAO),CNRFLO(AO),WNTF(AO),RPMNT(AO- 
1),DATOUN(9,60) ,FARN(A0),DANINP(15,60),TMTEMP(A0),TMPRES(A0),CORFLM- 
2(A0) ,DATOUM(9,60) ,DATMAC(A,60),DATALT(A,60),DE5LIM(15) 

COMMON /CENTER/ CGARMC60) 

COMMON /ARM/ ACCARM 

COMMON /DI5KK/ DI5KWI , ENGINE 

COMMON /SKIP/ SKIPIT 

^ DATA STORAGE DEFINITION ^ 

XXXXKXXXXKXXXXXXXXXXXXXXXXX 

LOGICAL PIKP, IPLT, 1510, ISI I, PLOT, SKIPIT 
INTEGER IDIDC60) , ILENG(AO) 

DIMENSION NUMNUMC17), IRNAMEC17), CORFLCCAO) 

NAMELIST / W/IWMEC,DESVAL,ACCS.IWT,IPLT,ISII,ISIO,IOUTCD,ILENG,DES- 
ILIM, I SCALE, SCALE, ACCARM, D I 5KWI, ENGINE, PLOT 

XKXXXXXXXXXXXXXXXXX 

K DATA STATEMENTS X 

XXXXXXXXXXXXKXKXXXX 

DATA IDUC,LSHAF,ENGU,SIU, IVALV/4HDUCT,AHSHAF,4HENGL,AHSIU ,4HVALV/ 


WEIGHT OF EACH COMPONENT 
ACTUAL LENGTH OF EACH COMPONENT 
ACCUMULATED LENGTH TO END OF COMPONENT 
RADIUS INNER INLET, OUTLET EACH STATION 
RADIUS OUTER INLET, OUTLET EACH STATION 
MECHANICAL DESIGN DATA OVERRIDES DEFAUL 
SHAFT DIAMETER INNER TO OUTER 
ACTUAL COMPONENT RPM 
WEIGHT ESTIMATION FLAG TRU= DO IT 
PLOTTER FLAG TRUE= DO IT 
GRAPHICS PLOTTING FLAG TRUE= DO IT 
ERROR FLAG 

OUTPUT UNITS O^ENGLISH, 0 SI 
INPUT UNITS 0=ENGLISH, 0 SI 

PRINT FLAG O^^SUMMARY, 1 =GENERAL , 2=DI AGNOSTIC 
COMPONENTS CONTRIBUTING TO OVERALL LENGTH 
FLAG = 0 COMPONENT NOT YET WEIGHED =1 YES 
BACKGROUND FLAG YE5=SKIP PROMPTS 


61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 <— 

72 

73 

74 

75 

76 

77 

78 <— 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 <<— 

108 

109 

110 
111 

112 <— 

113 

114 

115 

116 <— 

117 

118 

119 

120 
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DATA ILPC,IHPC,IFAN,IF0,IFI,IHPT,ILPT/3HLPC,3HHPC,3HFAN,2H^,2HFI 
13HHPT,3HLPT/ 

DATA PINP, YES/, TRUE. ,<*HY / 

TEST UTEST FLAG 

IF (IWT.EQ.O) GO TO 5<t0 
IF (IWAY.GE.O) GC TO 20 
DO 10 1=1, <i0 
WNTF(I)=WTF(I) 

FARN(I)=FAR(I) 

TNPRES(I)=TOPRES(I) 

TNTEMP(I)=TOTEMP(I) 

CNRFLO(I)=CORFLO(I) 

DAT0UN(1,I)=DAT0UT(1,I) 

DAT0UN(2,I)=DAT0UT(2,I> 

DAT0UN(8,I)=DAT0UT(8,I) 

DAT0UN(9,I)=DAT0UT(9,I) 

0 CONTINUE 

DO 15 1=1,900 

5 DANINP(I,1)=DATINP(I,1) 

0 DO 30 1=1,40 

CORFLCC I )=C0RFL0(I)/1. 5497255 
IF (JTYPE(I).EQ.l) II=I 

IF (TOPRESCD.GT.TMPRESCD) DATMACC 1 , 1 ) =DATINP ( 5 , II ) 

IF (TOPRES(I) .GT.TMPRESdn DATA L T( 1 , 1 ) =DATINP ( 9 , 1 1 ) 

IF (TOTEMPd) .GT.TMTEMPCin DATMACC 2 , 1) =DATINP ( 5 , 1 1 ) 

IF (TOTEMPd) .GT.TMTEMPd)) DAT ALT ( 2 , 11 =DATI NP ( 9 , 1 1 ) 

IF (DAT0UT(2,I).GT.DAT0UM(2,in DATMAC ( 3 , I ) =DATINP ( 5 , II) 

IF (DATOUT(2,I) .GT.DAT0UM(2,I)) DATALT ( 3 , I ) =DATINP( 9, II ) 

IF (CORFLCd).GT.CORFLMd)) DATMACC 4 , 1 ) =DATINP ( 5 , 1 1 ) 

IF (CORFLCCD.GT.CORFLMCD) DATALT C 4 , 1 ) =DATINP C 9 , 1 1 ) 

IF (TOPRESd) .GT.TMPRESCIl) TMPRES Cl ) =T0PRES Cl ) 

IF CTOTEMPCn.GT.TMTEMPCD) TMTEMPC I ) =TOTEMPCI ) 

IF CCORFLCCD.GT.CORFLMCD) CORFLMCI ) =C0RFLC( I ) 

IF CDAT0UTC2,!) .GT.DAT0UMC2,D) DAT0UMC2, 1)=DATOUT(2, 1) 

IF C JTYPECI) .NE.5) GO TO 30 
DATTRQ=DAT0UTC1 , I)/DAT0UTC2,I) 

IF (DATTRQ.GT.DATOUMCl,!)) DATOUMC 3 , 1 ) =DATOUTC 2, 1 ) 

IF CDATTRQ.GT.DATOUMCl,!)) DATOUMCl , I ) =DATTRQ 

THE FOLLOWING MAY BE USEFUL FOR DEBUGGING AND HAS BEEN LEFT IN 

WRlTEdO, 992)JTYPECI) 

992 FORMATC//' COMP TYPE'I3) 

WRITEC10,997) 

997 FORMATC DES TOTPRES N D TOTEMP’) 

URITEC10,999)TOPRESCI) ,TNPRESCI) ,TOTEMPd),TNTEMPd) 

WRITEC10,996) 

996 FORMATC* DES CORFLO N D M WTF’) 

WRITEC 10, 999 )CORFLOCI ) ,CNRFLOCI) , CORFLMCI ),WNTF( I) 

WR1TEC10,993) 

993 FORMATC DES EFF N D PR') 

WRITEC10,999)DATOUTC8,I) ,DAT0UNC8, I) ,DATOUTC9,I) ,DAT0UN(9,I) 
WRITEC10,995) 

995 FORMATC MAX TOTPRES N M TOTEMP*) 

WRITEC 10,999 ITOPRESCI) ,TMPRESC I) ,TOTEMPCI),TMTEMPCI) 

WRITEC10,994) 

994 FORMATC MAX RPM NMD*) 

WRITE(10,999)DATOUTC2,I) ,DAT0UMC2,I),DAT0UNC2,I) 

WRITEC10,991) 

991 FORMATC MTRQ N M TRPM MTRPM *) 

WRITECl 0, 999 IDATTRQ, DATOUMCl, l),DATOUT( 2,1), DATOUMC 3, I) 


121 

122 

125 <- 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 k 

176 

177 

178 

179 

180 
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C 999 

F0RMATC9F12.3) 

181 

C 

WRITECIO, 999)DATMAC(1/ I) ,DATNAC(2, I ) , DATMACC 3 , I ) 

132 

C 

WRITE(10,999)DATALT(1,I),DATALT(2,I),DATALT(3,I) 

183 

30 

CONTINUE 

189 


IF (IWT.GE.2) GO TO 90 

185 

C 

ZERO OUT OUTPUT ARRAYS 

186 


GO TO 590 

187 

90 

JSCALE=0 . 

188 


I5AVE=IWT 

189<— 


IF (IWT.NE.9) GO TO 70 

190 

50 

JSCALE=J5CALE+1 

191 


IF (J5CALE.GT.I5CALEC2)) GO TO 590 

192 


IF CJSCALE,GT.1.AND.ISCALE(1).EQ.2) IOUTCD=l 

193 


5CALEF=SCALE( J5CALE) 

199 


IF (JSCALE.GT.l) SCALEF=SCALE( J5CALE)/SCALE( JSCALE-1) 

195 

55 

IFCREVISE. EQ.YES) SCALEF=1.0 

196<— 


DO 60 1=1,90 

197 


WNTF(I)=WNTFC I )KSCAL EF 

198 


CMRFLO(I)=CNRFLO(I)^SCALEF 

199 


CORFLM(I)=CORFLM(I)KSCALEF 

200 


IDID(I)=0 

201 


DATOUNCl, I)=DATOUN(l, DXSCALEF 

202 

60 

CONTINUE 

203 


IF (JSCALE.GT.l. AND. IWT.EQ. 9) GO TO 120 

209 


IFCREVI5E.EQ.YES) GO TO 115 

205<~ 

70 

DO 80 1=1,5 

206 

80 

D5HAF(I)=0. 

207 


DO 90 1=1,60 

208 


WATE(I)=0 

209 


NSTAG(I)=0 

210 


IDID(I)=0 

211 


RPMT(I)=0. 

212 

90 

ALENG(I)=0 

213 


DO 100 1=1,90 

219 


TLENG(I)=0 

215 


RI(1,I)=0 

216 


ILENG(I)=0 

217 


RI(2, I)=0 

218 


RO(1,I)=0 

219 

100 

R0(2,l)=0 

220 

C 


221 

c 

NAMELIST READ OF WTEST DATA 

222 

c 


223 


CALL NAMEPR ( 9 , 1 0 , 8 , P INP ) 

229 


READ (8,W) 

225 


DO 110 1=1,60 

226 


DO 110 K=l,17 

227 


IF (IWMEC(1,I).EQ.NUMNUM(K)) IWMECC 1 , I ) =IRNAME(K) 

228 

*110 

CONTINUE 

229 


GO TO 120 

230<— 

115 

WRITE(20,720) 

231< — 


CALL NAMEPR(20,10,8,PINP) 

232<— 


READ(8,UJ) 

233< — 

C 


239 

c 

PROCESS LENGTH CONTRIBUTING VECTOR EXCEPT DUCTS AND SHAFTS 

235 

c 


236 

:i2o 

DO 290 1=1,90 

237 


NC=ILENGCI) 

238 

c 

WRITE(10,7777)NC,ILENG(I),JTYPE(NC> 

239 


IF (NC.EQ.O) GO TO 250 

290 



JT=JTYPECNC) • 2^1 

GO TO (2^0a60,200,l30,1^0,210,l90,150,180,2<i0,2<i0,2<*0,2'*0,2<>0),JT 2^2 

C COMPRESSOR 243 

130 CALL COMP (NO 244 

GO TO 230 245 

C TURBINE 246 

140 CALL TURB (NO 247. 

GO TO 230 248 

C MIXER 249 

150 CALL WMIXR (NO 250 

GO TO 230 251 

C PRIMARY BURNER 252 

160 IF (IUNEC(l,NO.EQ.IDUO GO TO 170 253 

IF (IWMEC(l,NO.EQ.IVALV) GO TO 220 254 

CALL COMBWT (NO 255 

GO TO 230 256 

C— DUCTS • 257 

170 CALL DUCTW (NO 258 

GO TO 230 259 

C NOZZLES 260 

180 CALL WTNOZ (NO 261 

GO TO 230 * 262 

C SPLITTER 263 

190 CALL WSPLT (NO 264 

GO TO 230 265 

C TRANSFER DIMENSIONS FOR WATER INJECTION 266 

200 CALL DUMMY (NO ' . 267 

GO TO 230' 268 

C HEAT EXCHANGER WEIGHT 269 

210 CALL HMEC (NO 270 

GO TO 230 271 

C VALVES 272 

220 CALL VALVWT (NO 273 

C ACCUME LENGTH 274 

230 IUP=JCONF(NC,l) 275 

IDN=JCONF(NC,3) 276 

TLENG(IDN)=TLENG(IUP)+ALENG(NC) 277 

IF (JT.EQ.6) TLENG(IDN)=TLENG(IUP) 278 

ID2=JC0NF(NC,4) 279 

IF (ID2.GT.0) TLENG(ID2)=TLENG(IUP)+ALENG(NC) 280 

IF (JT.EQ.6) TLENG(ID2)=TLENG(IUP) 281 

IU2=JC0NF(NC,2) 282 

IF (IU2.GT.0) TLENG(IU2.)=TLENG(IUP) 283 

IDID(NC)=I 284 

240 CONTINUE 285 

C LAST COMPONENT WAS A NOZZLE SET ENGINE MAXIMUM LENGTH 286 

250 ENGLEN=TLENG(IDN) 287 

C 2S8 

"c PROCESS REMAINING COMPONENTS 289 

C 290 

DO 340 1=1,60 291 

C WRITEC10,7777)NC,I,IDID(I) 292 

IF (IDID(I).EQ.l) GO TO 340 293 

C 294 

C— PROCESS COMPRESSORS, TURBINES, MIXERS, BURNERS, SPLITTERS 295 

. NC=JTYPE(I) 296 

. IF (NC.LE.O) GO TO 340 297 

GO TO (340, 290, 330, 260, 270, 320, 310, 280, 340, 340, 340, 340, 340. 340), NC 298 
C— COMPRESSOR 299 

260 CALL COMP (I) 300 
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GO TO 3<i0 i 301 

C-- TURBINES 302 

270 CALL TURB (I) 303 

GO TO 3A0 30A 

C— MIXER 305 

280 CALL UMIXR (I) 306 

GO TO 3^0 307 

C— BURNERS 308 

290 IF (lUMECd.D.EQ.IDUC) GO TO 3A0 309 

IF (IWMEC(l.I) .EQ.IVALV) GO TO 300 310 

CALL COMBWT (I) 311 

GO TO 390 312 

C VALVES 313 

300 CALL VALVWT (I) 319 

GO TO 390 315 

C~— SPLITTER 316 

310 CALL W5PLT (I) 317 

GO TO 390 318 

C HEAT EXCHANGERS 319 

320 CALL HMEC (I) 320 

GO TO 390 321 

C TRANSFER DIMENSIONS FOR WATER INJECTION 322 

330 CALL DUMMY (I) 323 

390 CONTINUE 329 

C URITE(10,7771)JTYPE 325 

C 326 

C PROCESS DUCTS 327 

DO 350 1=1,60 328 

IF nDID(I).EQ.l) GO TO 350 329 

NC=JTYPE(1) 330 

IF (NC.NE.2) GO TO 350 331 

IF (IWnEC(l,I)..NE.IDUC) GO TO 350 332 

CALL DUCTW (I) 333 

350 CONTINUE 339 

C U'RITE( 10,7771 )JTYPE 335 

C PROCESS NOZZLES ' 336 

DO 360 1=1,60 337 

IF (IDID(I).EQ.l) GO TO 360 338 

NC=JTYPE(I) 339 

IF (NC.NE.9) GO TO 360 390 

CALL WTNOZ (I) 391 

360 CONTINUE 392 

C WRITE(10,7771)JTYPE 393 

C ACCUME LENGTH 399 

DO 380 1=1,90 395 

C WRITE(10,7777)NC,I,IDID(1) 396 

IF ( JTYPE(I) .LE.O) GO TO 370 397 

IF (IDID(I).EQ.l) GO TO 380 398 

NC=JTYPE(I) 399 

GO TO (380, 370, 370, 370, 370, 370, 370, 370, 370, 380, 380, 380, 380, 380), NC 350 
370 IUP=JCONF(I, 1) 351 

1U2=JC0HF(I,2) 352 

IDN=JC0NF(I,3) 353 

ID2=JC0NF(I,9) 359 

TLENGC IDN)=TLENG( lUP ) + AL ENG( I ) 355 

IF (NC.EQ.6) TLENG(IDN)=TLENG(IUP) 356 

IF (IU2.GT.0) TLENG(IU2)=TLENG(IUP) 357 

IF (ID2.GT.0) TLENG(ID2)=TLENG(IDN) 358 

IDID(I)=1 359 

380 CONTINUE 360 
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i 


C PROCESS SHAFTS 361 

DO AOO J=l,5 362 

DO 390 1=1,25 363 

NC = KKINDS(ll,n 364 

IF (NC.LE.O) GO TO 400 365 

IF (IIJMEC(1,NC) .NE.L5HAF) GO TO 390 366 

IF (IWMEC(2,NC) , EQ. J) CALL SHAFT (NO 367 

90 CONTINUE 368 

00 CONTINUE 369 

370 

FIND ENGINE MAXIMUM RADIUS 371 

XR=0 372 

DO 420 I=1,N0STAT 373 

IF (XR.GE.ROd, D) GO TO 410 374 

XR=R0(1,I) 375 

10 IF (XR,GE.R0(2,D) GO TO 420 376 

XR=R0(2,I) 377 

20 CONTINUE 378 

379 

GET ENGINE TOTAL WEIGHT AND ALENG CONVERSION 380 

WATENG=0 381 

IF (ACCS.EQ.O) ACCS=.l 382 

WAT=0. 383 

DO 430 1=1,60 384 

IF (JTYPE(I).EQ.9) GO TO 430 385 

WAT=WATE(I)+UAT 386 

30 CONTINUE 387 

UATACC=ACCSXWAT 388 

IF (lOUTCD.GT.l) WRITE (10,580) lOUTCD 389 

IF (ISIO) WATACC=WATACC«C0NVER(3) 390 

IF (lOUTCD.GT.l) WRITE (10,610) WATACC 391 

DO 440 1=1,60 392 

WFACTR=1. 393 

IF (DESVAL(15,I) .NE.O. ) WFACTR=DESVAL(15, I) 394 

WATE(I)=WATE(I)XWFACTR 395 

IF (.NOT. ISIO) GO TO 440 396 

WATE(I)=WATE(I)XC0NVER(3) 397 

ALENG(I)=ALENG(I)XCONVER(l) 398 

40 WATENG=WATENG+WATE(I) 399 

400 

CONVERT RADIAL DIMENSIONS AND TLENO 401 

IF (.NOT. ISIO) GO TO 460 402 

DO 450 I=1,N05TAT 403 

RI(1,I)=RI(1,I)XC0NVER(1) 404 

R1(2,I)=RI(2,I)XC0NVER(1) 405 

RO(l,I)=RO(l,I)XCONVER(l) 406 

R0(2, I)=R0(2,I)XC0NVER(1) 407 

TLEN6(I)=TLENG(I)XC0NVER(1) 408 

50 CONTINUE 409 

410 

WRITE COMPONENT WEIGHT INFO 411 

460 UNITSI=ENGU 412 

IF (ISII) UNITSI=SIU 413 

UNITSO=ENGU 414 

IF (ISIO) UNITSO=SIU 415 

WRITE ( 10 , 590 ) UNITSI,UNITSO 416 

IF (IWT.EQ.4) WRITE ( 10 , 550 ) 417 

IF (IWT.EQ.4) WRITE ( 10 , 560 ) 418 

IF (IWT.EQ.4) WRITE ( 10 , 570 ) SCALE( JSCALE) 419 

IF (IWT.EQ.4) WRITE ( 10 , 560 ) 420 
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ORiGSNAL PAGE IS 
* OF POOR QUALITY 


IF nOUTCD.LT.l) GO TO 510 <>21 

WRITE nc, 600 ) ^22 

HLENG=0. ^23 

CGLENG=0. ^24 

CGWATE = 0. . <t25 

CGTOTM-0. <»26 

CGCOMP = 0. 627’ 

DO 500 1=1,60 A28 

NC=JTYPE(I) «29 

IF (NC.LE.O) GO TO 500 ^30 

C WRITE(10,7777)JT 631 

GO TO (670, 670, 670. 670, 670. 670, 670, 670, 670, 500, 670, 500, 500, 500), NC 632 

670 IUP1=JC0NF( 1,1) 633 

IUP2=JC0HF( I ,2) 636 

IDH1=JC0HF(I,3) 635 

IDH2 = JC0MF( 1 , 6) 636 

WRITE (10,620) I,WATE(I),ALENG(I),TLENG(IDN1),RI(1,IUP1),RO(1,IOP1- 637 

1) ,RI(1,IUP2),R0(1,IUP2),RI(2,IDN1),R0(2,IDN1),RI(2,IDN2),R0(2,IDN2- 638 

2) ,NSTAG(I) 639 

IF ( ALENG( I ) . EQ . 0 . ) GO TO 680 660 

HLEfiG = ALEHG( I )/2 . 661 

IF (IWMEC(l,I).EQ.ILPC.OR.IWMEC(l,l).EQ.IHPC) CG*X=ALENG( I )-CGARM( I- 662 

1) 663 

IF (IWMECa,!) .EQ.IFAN.OR.IWMEC(l,I).EQ.IFO) CGX = AL ENG( I ) -CGARMC I ) 666 

IF (IWMECd,!) .EQ.IFI.OR.IWNEC(l,I).EQ.IHPT) CGX = AL ENG( I ) -CGAR«( I ) 665 

IF ( IWMECC 1 , I ) . EQ . ILPT) CGX = AL ENG( I ) -CGARMC I ) 666 

IF ( IWMECC 1 , I ) . EQ . ILPC. OR. IWMECC 1 , I) . EQ. IHPC) HLENG=CGX -667 

IF (IWflEC-C i , I ) . EQ . IFAN.OR . IWMECCl , I) . EQ. IFO) HLENG=CGX 668 

IF (lUMEC(l.I) .EQ.IFI.OR.IWMEC(l,I).EQ.IHPT) HLENG=CGX 669 

IF ( IWMECC 1 , I) . EQ . ILPT) HLENG=CGX 650 

CGLENG=TLENG(IDN1 )-HLENG 651 

GO TO 6 9D 6'52 

680 IF (IWMECCl,!) .EQ.ISHAF) CGL ENG=CGARM( I ) 653 

IF ( IWMECC 1 , I) .NE. ISHAF) GO TO 500 656 

690 CGWATE=WATE(I) 655 

CGCOMP=CGWATExCGLENG 656 

C6T0TM=CGT0TM+CGC0MP 657 

500 CONTINUE 658 

CEHGRA=(CGT0TM+(WATACCXACCARM) )/WATENG 659 

C 660 

C MAKE SUMMARY PRINT 661 

510 IF (.NOT. 1310) GO TO 520 662 

XR=XRVC0NVER(1 ) 663 

ENGLEH=ENGLENHCONVER(l)' 666 

CENGRA=CENGRAXC0NVER( 1) 665 

520 WRITE (10,630) WATENG.WATACC, ENGLEN.XR 666 

IF CSCALECJSCALE) .EQ.l. ) SEXP01=WATENG 667 

SEXP02=1. 663 

IF (SCALE(JSCALE) .NE.l. .AND.IWT.EQ.6) SEXP0E=AL0GCWATENG/SEXP01 )/A- 669 

ILOGCSCALEC JSCALE)/SEXP02) 670 

IF CSCALECJSCALE) .EQ.l. ) SEXPOE=l, 671 

WRITE (10,660) CENGRA 672 

IF (IWT.EQ.6) WRITE (10,650) SEXPOE 673 

IF (JSCALE.GT.1.AND.ISCALE(1).EQ.2) GO TO 530 676 

IF (ENGINE. EQ. 2. ) GO TO 530 675 

IF (IPLT) CALL ENGPLT (ENGLEN.XR) 676 

IF(PLOT) GO TO 525 677< 

GO TO 530 678< 

525 IF (SKIPIT) CALL EGPLOT C ENGL EN , WATENG) 679< 

IF (SKIPIT) GO TO 530 68(K 



530 


5^i0 


C 

C 

550 

560 

570 

580 

590 

600 


610 

620 

630 


660 

650 

700 

710 

720 

750 


WRITE(20,700) 

READ(20,750) ANSWER 
IF(AN5WER.NE.YE5) GO TO 530 
CALL EGPLOT(ENGLEN,WATENG) 
WRITE(20,710) 

READ(20,750) REVISE 
IFCREVISE. EQ.YES) IWT=2 
ENGINE=1. 

IFCREVISE. EQ.YES) GO TO 55 
IWT=ISAVE 

IF (IWT.EQ.6) GO TO 50 

IWT = 0 

I0UTCD=2 

RETURN 


FORMAT (IH /26H ENGINE SCALING DATA ) 

FORMAT 

FORMAT (16H SCALE FACTOR ,F5.2) 

FORMAT (IH /16H xxxxxxxxxxxxx/lAH « ACCS WT 

1/16H ^ J</13H xxxxxxxxxxxx , II ) 

FORMAT (1H1,26H WEIGHT INPUT DATA IN ,A6,6H UNITS/27H WEIG- 

IMT OUTPUT DATA IN ,A6,6H UNITS//) 

FORMAT (69H COMP WT COMP ACCU UPSTREAM RADIUS DOWNS- 

ITREAM RADIUS /77H NO EST LEN LEN RI RO RI RO - 

2 RI RO RI RO NSTAGE/) 

FORMAT C/,11H ACCS WT=,F8.3) 

FORMAT (I7,F6.0,F7 .0,F6. 0,6 F5. 0,F6.0,3F5. 0,18) 

FORMAT C/,27H TOTAL DARE ENGINE WEIGHT=,F6.0,2X,12HACCESSORIES=,F- 
17.2,2X,23HESTIMATED TOTAL L ENGTH= , F6 . 0 , 2X , 25HEST IMATED MAXIMUM RAO~ 


2IUS=,F5.0) 
FORMAT (30H 
FORMAT (39H 
FORMAT (66H 
FORMAT (56H 
FORMAT (r2H 
1SVAL(1,6)=^^X 
FORMAT (A6) 
END 


ESTIMATED CENTER OF GRAVITY= , F6 . 0 ) 

ESTIMATED AIRFLOW SCALING EXPONENT IS,F6.3) 

DO YOU WISH A GRAPHICS PICTURE? YES=Y ; NO=N) 

DO YOU WISH TO MAKE CHANGES TO THE INPUT? YES=Y;NO=N) 
ENTER DESIRED CHANGES FROM TERMINAL , I . E . /39H 8W DE- 
^(DATA)XXXX &END) 
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SAMPLE CASE TERMINAL LISTING FOR INTERACTIVE MODE 


tu 

IWT = <( , IPLT = T,PLOT = T, 

ISIl^F, I5IO=F, IOUTCD=2 , 

UfHG(l) = 2.3,6,7,5,17,9,10.U,12,13.1<i, 
lWnEC(l,2):'FAN ' ,0,0, 9, 3X0, 

]HMtC(l,3)='SPLT’ ,6K0, 

1WMEC(1,9)='DUCT' ,2,5X0, 

IWMECa , 5) = 'NOZ ’,1,5X0, 
lWMEC(l,6)='tPC '1,5X0, 

H!HEC( 1 ,7)='DUCT ',2,5x0, 

1WMEC(1,8)=’HPC ',1,5X0, 

1U:1ECC1,9) = 'PBUR' ,1,5X0, 

IWMECC1,10) = 'HPT ' ,0,R,-S,3X0, 

ItJHECd , 11 > = 'DUCT' ,2, 5X0, 

IWMEC(1,12)=’LPT ',1,2, 6, 3X0, 

IUMEC(1,13)='DUCT' ,1,5X0, 

IWMECd, 19) = 'N02 ',1,5X0, 
lUHECd, 15) = ' SUAF ',1,12,3X0,6, 

IWHECd , 10)='SHAF' ,2,10,3X0,8, 

IWHF.Cd,17) = 'DUCT' ,2,5X0, 

DESVALd,2) = .55,1.8,.37,1.5,A.5,<f.5,.5,0.,0.,l.,0.,2.,l., 
DE5VAL(l,3)=15xo. , 

'•^DESVALd ,A) = .5,1 ■ , 0.,A., 11XQ., 

^ DE5VALd,5)=l.,lAx0., 

>DESVAtd,6) = .5,1.15 , .8, 1.5,2-. , 2 . , . A5 , 0 , 0 , 1 , 0 , 3 . , 1 . , 
DESVALd ,7) = .95,5.5nrT78. ,11X0., 

DE5VALC 1,8) =.95,1.31, .97,1.5,9.5,1.5,.3,0,0,1.,0.,3.,1., 
DESVALd, 9 ) = 100. , . 0 180 , 1 3X0 . , 

DESVALd,10)=.9,.20,1.5,1.5,1.5,.5,150000.,3.,l.,6xO., 
DESUALd,ll) = .5,3.5,0.,-l.,llxO., 

DESVALd,12)=.5,.327,1.5,2.,9.,.6,150000.,l.,l.,6XO., 
DE5VALd,13) = .6,.5,0.,19.,llX0., 

DESVALd,19)=l.,19X0. , 

DESVALd . 15) =50000 . ..3, .85, 12X0., 

Dn5VALd, 16)=50 000. , .3,13X0 ., 

DESVALd , 17)= .3,8. 9, 0 11x0., 

SEND 


AS0iWATE2IT TXXXX, RET=T , PRINT=NO , NEWLIB= JOB . BODLIB 

DEFAULTS ARE IH=WHI TLOW , OUT=ANSWER . IN , PRINT=YES , MAPS=ADV3DMAP , L IMEN=W, NEWL IB=NONE, RET=P , BACK=NO 
THIS RUN, IH=CF6 ,OUT=ANSUER.CF6,PRINT=NO,MAPS=KAPS.PWSCAR,LIMEN=W,NEWLIB=JOB.BOBLIB,RET=T,BACK=HO 
JOB.BODLIB LOADED 
BLOCK DATA BLOCKWT IS LOADED 
ENGWT LOADED 
00001 
00002 
00003 
00009 
00005 
000D6 

00007 

00008 

LAST QUALIFIER IS FROM NEPCAL 

00009 

00010 
0001 1 
00012 
00013 
00009 

EXECUTION HAS BEGUN 
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00005 

> > > > > SUCCESSFUL CONVERGENCE < << < < 

N1T = 0 

DO YOU U15H A GRAPUICS PICTURE! YES=Y ; KO=M 

Y 

DO YOU WISH MULTIPLE FRAMES? YES=Y ; NO=N 
N 

DO YOU WANT COMPONENT OUTPUT? YES=Y ; NO=N 

Y 

NOT ENOUGH ROOM BETWEEN LABELS. NOLAB REDUCED BY 1/2. 

GRAPHICS DEVICE NOT DEFINED BY DDEF. 

ENTER UNIT NAME, DEFAULT TO CANCEL. 

LAOOl 

MAKE SURE GRAPHICS DEVICE STILL CONNECTED TO TEMPO/5035 
DO YOU WISH TO MAKE CHANGES TO THE INPUT? YES=Y;NO=N 

Y 

00012 

WATENG=+.E7A8691AE40A 
WATACC-+.8205A126E+05 
WATE - 

(1) + . OOOOOOOOE+00 + , 20693882E+0A + , 00 0 0 0 0 0 0E+ 0 0 + .21138219E+02 

(8) 4.10A29136E+0A + . A6836670E+03 + . 278<i5117E + 03 +. 1<»6098A9E+02 

(15) + . A6A37256E405 + . 65A9 I 1 96 E+02 + .59371A90E+02 + . OCOOIOOOE+00 

(22)'*(60) CONTAINS 4 . OOOOOOCOE+00 
ENTER DESIRED CHANGES FROM TERMINAL , I . E. 

l\4 DESVAL(l,6) = >f><^>?(DATA)X)<K<< &END 
8U DESVAL(3,6)=.8A,DE5VAL(3,A)=35. »0, iEND 
DO YOU WISH A GRAPHICS PICTURE? YES=Y ; NO=K 

Y 

DO YOU WISH MULTIPLE FRAMES? YES=Y ; NO=N 
N 

DO YOU WANT COMPONENT OUTPUT? YES“Y ; HO=H 

Y 

NOT ENOUGH ROOM BETWEEN LABELS. NOLAB REDUCED BY 1/2. 

DO YOU WISH TO MAKE CHANGES TO THE INPUT? YES=Y;NO=N 
N 

00012 

WATENG=4.85592891E40A 
WATACC-4,7950A028E403 
WATE = 

(1) 4.GOCOOOOOE400 4.20693882E40A 4 . OOOOOOOOE400 4 . 21235367E4Q2 

(8) 4.10929136E40<< 4 . <i6S36670E403 4 , 278<* 51 17 E403 4 . I^609899E402 

(15) 4.<*5298218E403 4 . 65991 1 96 E4Q2 4 . 3937 19 90 E4 02 + . OOOOOOOOE400 

(22):(60) CONTAINS . 0 00 0 0 00 0E40 0 
DO YOU WISH A GRAPHICS PICTURE? YE5=Y ; NO=N 
N 

00012 

WATENG=4,65969336E409 
WATACC=4.60576050E403 
WATE = 

(1) 4.00000000E400 4.16723853E409 4 . OOOOOOOOE400 4 . 16509925E402 

(8) 4,79099707E403 4 . 38088965E+03 4 .20086632E403 4 , 1 0955938E402 

(15) 4.38575269E403 4 , 59026550E402 + .28171722E402 4 . OOOOOOOOE400 

(22):(60) CONTAINS + . o 0 0 0 C 0 00 E+0 0 
DO YOU WISH A GRAPHICS PICTURE? YES’Y ; NO=M 
N 

00012 

WATENG=4.)0718051E405 
WATACC=4.10039109E409 
WATE = 

(1) 4.00000000E+00 4 .2500R088EH 09 4 . O0OO0000E400 4 . 26188171 E+ 02 

(8) 4.U893208E409 4 . 55517700E403 +. 36383301 E4C3 4 . 1 9205093E402 

(15) 4.51733716E+03 4 . 7 191S152E^ 02 + . 51759623E402 4 . 0 0 0 0 00 00 £4 00 

(22):(60) CONTAINS + , OOOOOOOOE400 

00006 

PERPF( 9) =4. 51 90356 3923376270405 
PERPF( 5) =4.38855219022369210+00 
PERPF( 9) =4.51903563923376270405 
PERPFC 10) =4.38855219022369210400 
TERMINATED: EXIT IN USER PROGRAM 
CANCELLED: DUM.T2 UNKNOWN. 

A50: 


4.35629565E403 

4.29S69705E409 

4.00000000E400 


4.35957813E403 

4.19867976E409 

+.OOOOOOOOE400 


4.33938579E403 

4.1A909509E409 

4.00000000E400 


4.37977100E403 

4.25965910E409 

4.00000000E400 


4.12359275E409 

4.00000000E+00 

4.00000000E400 


4.U909098E409 

4.0000000CE+00 

4.00000000E400 


4-.10575210E409 

+.OOOOOOOOE400 

+.00000000E40O 


4.19757178E40^ 

4.00000000E400 

4.00000000E400 


4.18^32083E4C 

4.18697968E40 

4.00000000E40 


4.20906189E40 

4,29930663E40 

4-.OOOOOOOOE40 


4.19601971E40: 

4.19999562E40: 

4.00000000E40< 


4.26829696E402 

4.29916729E40: 

+.00000000E40C 
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APPENDIX E 


IBM-370 GRAPHICS SUBROUTINES (INCLUDES ONLY THOSE USED IN RUNNING THE 
PROGRAM) 

CHAUS — Print Character Data on a Plot 

The CHARS subroutine allows the user to print character data 
anywhere on a plot. The routine can be used to label a 
point. 


IGeneral Form | 

, 

(CALL CHARS (nchar, char, orient, X, y,isiz) I 


nchaf 

is an integer constant or variable the absolute value 
of which specifies the number of characters to be 
printed. For nchar < 0, the x- and y-coordinates of 
the starting position for the character data will be in 
user data units. For nchar > 0, the coordinates of the 
starting position for the character data will be in 
relative (0-10), units. 

char 

is the array of characters to be printed. The 
dimension of char must be large enough for the number 
of characters desired. The characters must be packed 
four to a word (A4 format) . See Appendix A for EBCDIC 
and integer equivalents for characters available on 
each device when a stroke table is not in use. 

orient 

is the angle of orientation in degrees from the 
horizontal. (OANGL from GRCOH will be added to this 
angle) . If the angle is 0 degrees, and OANGL is 0 
degrees, the characters will be printed right side up 
and horizontally. 


X 

is a floating point constant or variable which 
represents the x-coordinate of the starting position 
for the character data. If x is expressed in user 
units, NCHAR must be less than 0. If x is expressed in 
relative (0-10), units, NCHAR must be greater than 0- 

I 

is a floating point constant or variable which 
represents the y-coordinate of the starting position 
for the character data. If y is expressed in user 
units, NCHAF. must be less than 0. If y is expressed in 
relative units, NCHAR must be greater than 0. 
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isiz 

is an irit^qer constant or variable the value of which 
specifier, the size* of the characters to be printed. 
Its value must be in the range -127 to +127 and is 
interpreted as follows: 

< 0 Display italicized characters. 

Each character is drawn in a square whose size in 
internal units is calculated by the following formula: 

. box size = 128/(143-IABS <ISIZ) ) ♦ (TABS (ISIZ) +1) ♦15.0. 
To find the number of characters that will fit across a 
frame, divide the box size (found above) into 16384 or 
see Appendix B. 

Note: The standard character size is 15. 

Programming Note: If x and y are specified in user’s units, 

the user must have defined XMAX, YBAX, XBIN, and YMIH by 
previous calls to SCALE, AXIS or NAXIS.. If nonstandard axes 
are to be used and x and y are specified in user’s units, 
the user must have defined the plot boundaries by previous 
calls to AXIS or NAXIS. 

Example: The following is an example of using the CHARS 

subroutine to label a point. 

10 DIBENSIOK CK(3) 

20 DATA CH/» (3.0*, * ,5.0*,*) */ 

30 CALL CKAES(9,CH,0.0,3.0,U.0,15) 

10 Defines the dimension of the array which 

contains the characters to be plotted. 

20 Initializes the array with the information to be 

printed. 

30 Causes the printing of the label (3.0, 5.0) in 

standard size characters, upright and parallel 
to the X axis beginning at relative coordinate 
location (3. 0,4.0). 
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CORNER - Delete / Restore Corner Harks | 

The CORNER subroutine will remove or restore the corner | 
marks that define a frame. I 


{General Porm i 

, , 

{CALL CORNER (iopt) I 


iopt 

is an integer indicating whether to remove or restore 
the corner marks. 

Specified as: 1 - remove the corner marks 
0 - restore the corner marks 

Functional Description: The option will remain set until 

changed. This call should be made before the frame is 
displayed (i.e., before the applicable DISPLA call). 


DISPLA — Display 


The DISPLA subroutine defines the end of a display or plot 
and initiates transmission of orders to the deyice- 


IGeneral Form | 

I — 1 

ICALL DISPLA (OPT) | 


opt 

is a fullword integer constant or variable the value of 
vhich indicates the status of the buffers after the 
call. 

specified as either 0 or 1. 

OPT=0 specifies that the general orders buffers are not 
to be cleared following the display. 

OPT=1 specifies that the general orders buffers and 
GRCOM values will be cleared following the display* 
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GARC — Approximate circles and arcs. 

The GARC subroutine allows the user to approximate circles 
and arcs by points or vectors. 


IGeneral Form | 

, 1 

ICALL GARC (centX/Centy, rad, angle, iopt) | 


centx 

is a real constant or variable the value of which is 
the x-coordinate of the center of the circle/arc in 
relative units. 

• 

centy 

is a real constant or variable the value of which is 
the y-coordinate of the center of the circlt/arc in 
relative units. 

rad 

is a real constant or variable the valc» of which is 
the radius of the circle/arc in relative units. 

angle 

is a real array for specifying variable data to the ARC 
subroutine. Its entries should contain the following 
information: 

angle (1) The incremental angle for points or 
vectors. If angle (1) is 0, the default 
for the device is used. 

angle (2) The beginning angle for the arc in 
degrees. 

angle (3) The counterclockwise sweep angle for the 
arc. 

iopt 

is an integer constant or variable which indicates 
whether points or vectors are to be used to approximate 
the arc. Its value must be one of the following: 

0 - use default for device, 

1 - approximate with points. 

2 - approximate with vectors. 
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GPLOT 


Generate a Plot 


The GPLOT subroutine enables the user to plot single or 
multiple curves or lines as vector, point, symbol, or 
vector-symbol plots. 


(General Form I 

(CALL GPLOT (x,y,ivars) | 


X 

is a floating point array of x-coordinates of the 
points to be plotted or the array of radius-values for 
a polar plot. 

y 

is a floating point array of y-coordinates of the 
points to be plotted or the array of theta-values (in 
degrees) for a polar plot. 


ivars 

is an integer array for communicating other options to 
the subroutine. 


When plotting a single curve, the information supplied 
in ivars should be specified as follows: 


ivars (1) 
ivars (2) 
ivars (3) 


ivars (4) 
ivars (5) 

ivars (6) 

ivars (7) 


The number of elements in ivars. 

The number of points in the curve. 

The type of plot, same as in IVABS(2) 
for multiple curve plot. 

Default: 0 is assumed. 

The symbol code for a symbol plot. 

The symbol frequency for a symbol plot. 
Default: 1 is assumed. 

The size of the symbol to be plotted. 
Default: 15 is assumed. 

The interval option. 

0 - scaled (adjusted) to include MIN 

and MAX. 

1 - exact HIM and MAX interval. 
Default is 0 (adjusted) . 


Notes: If a vector plot is desired, ivars(3) through 
ivars (7) need not be specified. If a point plot is 
desired, ivars (4) through ivars (7) need not be 
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specified. 

For plotting multiple curves, the entries in ivars 
should contain the following information: 

ivars (1) The number of elements in ivars. 
ivars(2) The type of plot. Its value must be one 
of the following: 

0 vector plot 

1 point plot 

2 vector-symbol plot 

3 symbol plot 

4 polar vector plot with no call 
to PAXIS 

5 polar point plot with no call 
to PAXIS 

6 polar vector-symbol plot with 
no call to PAXIS 

7 polar symbol plot with no call 
to PAXIS 

12 polar vector plot with PAXIS call 

13 polar point plot with PAXIS call 

14 polar vector-symbol plot with 
PAXIS call 

15 polar symbol plot with PAXIS call 

16 log vector plot 

17 log point plot 

18 log vector-symbol plot 

19 log symbol plot 

20 log-polar vector plot with NAXIS 
call 

21 log-polar point plot with NAXIS 
call 

22 log-polar vector-symbol plot with 
NAXIS call 

23 log-polar symbol plot with NAXIS 
call 

28 log-polar vector plot with axes 

29 log-polar point plot with axes 

30 log-polar vector- symbol plot with 
axes 

31 log-polar symbol plot with axes 
ivars (3) The symbol code for a symbol plot (See 

Appendix A) . 

ivars (4) The frequency of printing the symbol for 
a symbol plot. 

Example: ivars (4) = 3 would cause the 

printing of the special symbol at every 
third point. 

Ivars (5) Specifies whether or not duplication of 
coordinates is to take place when 
plotting multiple curves. Its value 
must be one of the following: 

♦1 duplicate x-coordinates 
0 no duplication of coordinates 
“1 duplicate y-coordinates 
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ivars(6) The size o: the symbol to be plotted. 

Default: 15 is assumed. 
ivars(7) The interval option. 

0 - scaled (adjusted) to include MIN 

and MAX. 

1 - exact MIN and MAX interval. 

If not specified* an adjusted interval 
vill be created. 
ivars(8) - ivars (i vars (1) ) 

One entry for each curve specifying the 
number of points in that curve to be 
plotted. 


If there have been no calls to AXIS before PLOT- is called, 
the standard size axes with 90 degree orientation will be 
drawn with no titles, but with 10 intervals and 11 grid 
lines or tick marks and 11 labels, (tick marks will be the 
default) . 

TITLE can be called to print titles for standard x or y axes 
or to print a plot title. TITLE should be called following | 
GPLOT. I 

If XHAX and XHIN or YMAX and YMIN are equal, GPLOT assumes | 
maximum and minimum have not been determined and calls SCALE | 
to scan the array and determine the maximum and minimum 
coordinates. Depending on the interval option, if scaled or 
adjusted was requested, GINTVL will be called to create an 
interval containing the data MIN and MAX. 

To preset XHIN, XMAX, YMIN, YMAX, use the AXIS routines, 
SCALE, SCLBAK, or GINTVL. | 

Specifying a symbol type plot allows the user to plot 
symbols or characters at random locations. The symbol or I 
character will be centered at the X,Y location. j 

Note: In FOET5AN programming, selective array locations | 

cannot be defaulted. For example, if only array location 1, j 
2, and 3 in IVALS are specified, the remaining locations | 
will be defaulted. However, if location 7 is specified, | 
locations 1-6 must also be specified. Zero is not j 
recognized as a default value, j 

Example: The following example shows the use of the GPLOT 

subroutine to plot multiple curves as point plots. 

5 DIMENSION IVAHS(9) 

10 DIMENSION X(30),Y(60) 

20 XSTABT =0.0 

30 DO 40 1=1,30 

X(I) = XSTART 

Y(I) = XSTAET*XSTABT * XSTAET ♦ 1 
Y(I+30) = XSTART*XSTAPT - XSTAET ♦ 2 
XSTAET = XSTAET +0.5 
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40 CONTIHUE 

45 IVAP.5{i) = 9 

50 IVARS(2) = 1 

55 IVARS(3) = 0 

60 IVARS{4) = 0 

65 1VARS(5) = +1 

70 IVARS(6) = 0 

80 IVARS<7) = 0 

90 1VARS(8) = 30 

95 IVARS(9) = 30 

100 CALL GPLOT(X,y,IVARS) 

5 Dimensions the IVARS array. 

10 Defines the x and y arrays. There will be 2 

y~coordina tes for each x-coordinate. 

20 Initializes XSTABT at 0.0. 

30-40 This loop places values in the X and Y arrays. 

45 Defines the number of elements in IVARS to be 

seven. 

50 Specifies a point plot. 

55-60 Specifies no symbol code or index. 

65 Specifies that there will be duplicate 

2 -coordinates. 

70 Specifies no symbol size (no symbol) . 

80 Specifies interval option adjusted. 

90-95 Specify that there will be two curves with 30 
points in each. 


,OR5G'NAL PAGE IS 
OF POOR O'-'ALITY 
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GBIMIT — Graphics package initialization. 

The GRINIT subroutine initializes the graphics package. 


{General Porn 
{CALL GBINIT 



There are no parameters for this call. 

Functional description: GRINIT causes initialization of the 
graphics package and prompts the user for a device (if one 
has not been previously defined). This call is optional in 
that a call to any graphics subroutine will cause 
initialization if not previously accomplished. 
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NUMBEB — Convert an Integer or Real Number to a Character 
Array 

Tne NOMBEB subroutine allows the user to convert an integer 
or real number to an array of printable characters. 


I General Fora | 

jcALL NUMBER (ityp,rnum, nochar, nodec, char) j 


ityp 

is an integer constant or variable the value of which 
specifies the type of number to be converted. Its value 
oust be one of the following: 

0 Specifies an integer halfword is to be 

converted. 

1 Specifies an integer fullword is to be 

converted. 

2 Specifies a single precision real number is to 
be converted to •£* format. 

3 Specifies a double precision real number is to 
be converted to *D» format, 

• 4 specifies a single precision real number is to 
be converted to *F* format. 

rnum 

is a constant or variable (either integer or real as 
specified in ityp) the value of which is the number to 
be converted. 

nochar 

is an integer constant or variable the value of which 
specifies the number of characters desired in the 
output array. 

nodec 

is an integer constant or variable the value of which 
specifies the number of positions right of the decimal 
point. 

char 

is the output array of printable characters. Char must 
be specified large enough to hold the number of 
characters requested (nochar). The characters are 
packed four to a word (A4 format). 
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Note: If ITYP = 2 or 3, and NOCHAR is less than 8 or NODEC 

is greater than (NOCUAR-6) , or if ITYP =0, 1 or 4, and 

NOCHAR is not large enough to hold the number of significant 
digits in RNUH, or if ITYP is invalid, asterisks are placed 
into the output array. 

Example: The user wants to convert the real number XMAX to 
an array of printable characters in an array named CHARS. 
The printed number will be 8 digits long with 2 decimal 
places and will be in S format. 

.10 DIMENSION CHARS (2) 

20 CALL NUMBER (2, XMAX, 8, 2, CHARS) 

10 CHARS contains two elements (8 characters in A4 

format) . 

20 Calls the NUMBER subroutine: 

2 - convert a real number 

XMAX - contains the number to be converted 
8 - characters in output array 
2 - positions right of the decimal point 
CHARS - stores the printable characters 
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SETFBM — Set number of frames for plot 


The SETFBH subroutine allows 
of frames to be plotted. 

the user to 

specify the number 

{General Form 


1 

1 

{CALL SETFRM (frames^ 


j 

1 

frames 

is an integer constant 

or variable 

the value of which 


specifies the number of frames in the user's plot, 
creating a larger physical plotting area in the x 
direction. Its value may range from 1 to 21 inclusive. 

Example: The following is an example of using the SETFKM 
subroutine to set the number of frames to 2. 

30 CALL SETFEM C2) 

30 sets plot boundaries for a plot of 2 frames (up 
to 20. 0 relative units in x direction) . 
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TITLE — Print Plot or Axis Title 

The TITLE subroutine allows the user to print a plot title 
for a plot or to print titles for the axes used in the plot. 


IGeneral’ Form | 

, 

ICALL TITLE (ityp ,nchar,isiz, char) |- 


ityp 

is an integer constant or variable the value of which 
specifies what type of title is to be generated. Its 
value must be one of the following: 

■ 0 A plot title is to be printed, centered on the 
plot. 

1 A plot title is to be printed, left justified on 
the plot. 

2 A plot title is to be printed, right justified on 
the plot. 

3 A .y-axis title is to be printed- Mot valid for 
polar plots. 

ft An x-axis title is to be printed. Not valid for 
polar plots. 

nchar 

is an integer constant or variable the value of which 
specifies the number of characters in the title. 

isiz 

is an integer constant or variable the value of which 
specifies the size of the characters to be printed. 
Its value must be in the range -127 to +127 and is 
interpreted as follows: 

< 0 Display italicized characters. 

Bach character is drawn in a square whose size in 
internal units is calculated by the following formula: 
box size = (128/(143-IABS (ISIZ)) ♦ (lABS (ISIZ) +1) *15.0, 
To find the number of characters that will fit across a 
frame, divide the box size (found above) into 1638ft or 
see' Appendix B. - * 

Bote: The standard character size is 15. 
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char 

is an array of characters for the title. The dimension 
of CHAR must be la.rge enough for the number of 
characters desired. The characters oust be packed four 
to a word (A4 fornat). 

Programming Note: If TITLE is called to print plot or axis 

titles for a plot with nonstandard placement or orientation 
of axes, the plot boundaries must have been defined by a 
previous call to AXIS or NAXIS. 

Example: The following example shows hov the TITLE 

subroutine can be used to print axes titles and a plot 
title. 

10 DIMENSION XTITLE(6), YTITLE(6)., PTITLE (5) 

20 DATA XTITLE/'SAMP* , 'LE X*,*-AXI»,'S TI',‘TLE V 

21 DATA YTITLE/’ SAMP' , *LE Y','-AXI','S TI','TLE '/ 

22 DATA PTITLE/'SAMP' , 'LE P','LOT *,»TITL','E V 

30 CALL TITLE <4 , 19, 1 5, XTITLE) . 

HO - CALL TITLE (3 , 19 , 1 5 , YTITLE) 

50 . CALL TITLE (0, 17,35, PTITLE) 

10 Defines the arrays from vhich the titles will be 

printed. 

20 Initializes the arrays defined in 10 with the 

actual titles to be printed. 

30 Causes the X'-axis title to be printed in standard 
size characters. 

HO Causes the y-axis title to be printed in standard 
size characters. 

50 Causes the plot title to be printed ’ in large, 

(3/6"), characters and centered. 
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ORIGINAL PAGE IS 
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XAaIS, YAXIS, GAXIS — Draw Nonstandard Axes 


The axis subroutines can bo 
lengths, in nonstandard 
grid/tick nark frequency, 
degrees and SO degrees. It 

a range of plotted coordinate values. If no 
drawn the user must call NAXIS before PLOT is 


used to draw axes of nonstardard 
positions, with nonstan'lard 
and/or at angles other thatn 0 
also permits the user to specify 

axes are to be 
called. 


XAXIS is used to draw x axes, YAXIS is used to draw y axes, 
and GAXIS is used to draw additional axes. Both XAXIS and 
YAXIS should be called in succession to avoid defaulting the 
other's values. GAXIS will affect the scaling and 
placement of the user's data. 


(General Form 

, 

(CALL XAXIS (X, y, vars) 
jCALL YAXIS (x,y, vars) 
(CALL GAXIS (X, y, vars) 


is a real constant or variable which represents, for 
GAXIS calls, the x-coordinate of the starting position 
of the axis in relative units, and for XAXIS or YAXIS 
calls, the x-coordinate of the lower left plot 

boundary. A value of -1.0 represents the standard 
starting x-coordina-te position* 

NOTE 1; The values entered on the first call, (to 
either X or YAXIS) , should be the same as the values 

used for the second call, (to either X or YAXIS), If 

. they are not, the results will be unpredictable. 

NOTE 2 : If the standard X and Y location is specified 

(-1.0), the package will calculate the X and Y 

location, allowing room for horizontal y-axis labels. 
If the X and Y location is specified, it will be used 
even if there is insufficient room for the labels. 

y 

Same definition as for x above, except that this value 
represents a y-coordinate, 

vars 

is a real array for communicating other variables and 
options to the subroutine. The entries in vars contain 
the following information: 

vars(1) The number of elements in vars. (0-9) 
vars (2) The axis length in relative units, (0-10 
for Y; 0-210 for X). Default or -1. 
gives standard length (XUR-XLL) if XAXIS 
is called; (YUP.-YLL) if YAXIS is called. 


45 



VARS (3) 


▼ars (4) 


vars (5) 


vars (6) 


vars (7) 


vars (7) 


vars (8) 


Thp angle of orientation in degrees. 
Specify as 0 degrees for x-axis and 90 
degrees for y-axis. The default for 
vars (3) is 0 degrees + OANGL (from GKCOM) 
if xaxis is called; 90 degrees + OANGL if 
YAXIS is called. 

The value represented by starting point 
of axis (maximum or minimum). The 

default for vars (4) is (XMIN/YMIN) . | 

Note: This value must be entered unless | 
the user has already specified XMIN/YMIN | 
by using SCALE, SCLEAK, or GINTVL. The I 
value entered for vars (4) replaces the 

current XHIN/YMIN value in GRCOM. 

The value represented by end- point of 
axis (maximum or minimum) . The default 
for vars(5) . is (XMAX,YMAX). Note: The 

same rule applies as in note above for 
vars (4), except this value replaces 
XHAXIYHAX. 

The number of grid line/tick mark 
intervals. Its value can be one of the 
following: 

n the number of grid line/tick 
Bar); intervals desired. 

.5 standard number of grid 
line/tick mark intervals. (10 
intervals-^1 1 tick marks). This 
is the default value if not 
entered. 

-n the number of cycles for a 
logarithmic scale. 

If 0 is input, 1 interval with 2 grid 
lines/tick marks will be created and the 
default no of labels will be set to 0. 

If 1 is input, 1 interval with 2 grid 
lines/tick marks will be created and the 
default no of labels will be set to 2. 

Grid line/tick mark option (for XAXIS or 
YAXIS calls only). Its value should be 
one of the following: 

0 default for device, (this is 
the default if not input or if 
incorrect) , 

1 grid lines. 

2 tick marks. 

The associated YAXIS angle (for GAXIS 
calls only) .. Default is vars (3) ♦ 90 

degrees * OANGL. 

The number of labels per axis. If this 
element is not specified, or is negative, 
the number of labels will correspond to 
the number of grid lines/tick marks, 
(intervals ♦ 1), except when vars (6) was 
input as 0., in which case, the number of 
labels will be 0. The y-axis labels will 
be horizontal and on the counterclockwise 
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si do of th».' axis. 

vars(9) TJiO user data value (DX|Dlf) for 
displacement of this axis (on. the other 
axis) in user units. Default if not 
entered or invalid is zero. The DX|DY 
values will be calculated automatically 
for the user, (for XAXIS 0 YAXIS calls 
only), to cross the axes at 0,0 if -255. 
is entered. If this value is entered as 
-255 for a call to either X or YAXIS, it 
oust also be entered as -255 for the 
other call. 

wars (10) The size of the labels generated by this 
routine. Value from 0 to 127. Default if 
not specified will be 10. 

Note: In FORTEAN programming, selective array locations | 

cannot be defaulted. For example, if only array locations j 
1,2 and 3 are specified, 4-10 will be defaulted. However, j 
if 1 and 2 are specified and 6 and 7 are specified, 3, 4, j 
and 5 must also be specified (although 8, 9, and 10 can be j 

defaulted). Zero is not recognized as a default value. | 

Example: The following is an example of using the AXIS 

subroutines to draw axes with nonstandard specifications. 

The user wants the x axis to be drawn at a 0 degree angle 
with the horizontal, the y axis to be drawn at a 90 degree 
angle with the horizontal, and the origin to be at user data 
point 0.,0.. 

DIBENSIOK VARS (9) 

10 VARS(1) = 9. . 

20 VARS (2) = 7. 

30 VARS (3) =0. 

40 VARS (4) = -20. 

50 VARS (5) =20. 

60 VARS (6) = 4. 

70 VARS (7) =2. 

80 wars (8) = - 1. 

90 wars (9) = -255. 

100 CALL XAXIS (1., 1., VARS) 

110 VARS (3) =90. 

120 VARS (4) = -100. 

130 VARS (5) = 200. 

140 vars(6) =6. 

150 CALL YAXIS (1., 1., VARS) 

10 Defines the number of elements in VARS to be 9. 

20 Sets the length of the x axis equal to 7 

relative units, and x upper right plot boundary 
to X ♦ length or 8 relative units. 

30 Sets the angle of orientation for the x axis to 
0 degrees. 

40 Sets the value represented by the starting point 
of the X axis at -20. 

50 Sets the value represented by the ending point 
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GTE5H Command 


This command is required to close the graphics data set and 
perform any necessary cleanup functions. 


I Operation I Operand | 

, , , 

IGTERM 1 I 


Note: there are no operands. 

Functional Description: GTERM initiates the closing of the 
graphics data set, and any necessary cleanup functions. For 
offline devices, GT2FH establishes the BWQ task to transfer 
the dr.taset to tape for plotting. GTERM also frees virtual 
memory in which graphics entities have been saved. ’ 
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PROCEDURE DEFIMITION: WATE2IT 


ORIQ5NAL PAGE IS . 

OF POOR QUALITY 

PR0CLI5T WATE2IT 
WATE2IT OOOQOOO PROCDEP WATE2IT 

WATE 2 IT 0000100 PAR AM IN=$ IH , OUT=$OUT , PRINT=$PRINT ,MAP 5 =$riAPS , L IMEN=$L IMEN , NEWL IB=$NEWt IB , RET=$RET ,BACK=$BACK 
UATE2IT 0000200 D ’DEFAULTS ARE IN=1JHITIOW»OUT=AKSWER . IN, PR1NT=YES»MAP5=ADV3DMAP, LIMEN=W,HEWLIB=NONE, RET=P,BACK=HO* 
WATE2IT 0000300 IF ’ $OUT * -= ’ NONE’ ; “ 

WATE21T OOOOAOO D ’THIS RUN . , OUT=ANSWER . $0UT , PRINT=$PRINT ,«APS=0MAPS , L IMEN=SLII1EN, NEWLIB=$NEWL IB ,RET=$RET ,BACK=$BA{ 
WATl2!T 0000500 IF ’ $OUT * = ’ HONE ’ ; - 

WATE2I 1 0000600 D ’THIS RUN , IN = SIN . 0UT = AN5UER . BIN , PRINT =$PRINT . MAPS=$MAPS , I IMEN=$LIMEN, NEULIB=$NEWLIB ,RET=$RET . BACK=BBACJ 

WATE2H 0000700 DEFAULT LIMEN^BLIMEN 

WATE2IT 0000800 IF ’ $OUT NONE ERASE ANSWER. $OUT 

U'AT£2IT 000090C IF * SOUT* = ’NONE’ ; ERASE ANSWER. BIN 

WATEZiT OGOIOOO RELEASE FTO 1 FO 0 1 ; DDEF FT 0 1 FO 0 1 , VS , DUM . T1 , RET = T 

WATE2I1 OOOIlOO RELEASE FT02F00 1 ; DDEF FT02FC 0 1 , VS , DUM . T2 , RET=T 

tM7E3IT 000120C RELEASE FT03E001 ;DDEF FT03F00 1 , VS , DUM . T3 , RET-T 

WATE2IT 0001300 RELEASE FTOAFOOl ; DDEF FT 0 AFOO 1 , VS , DUM . TA , RET = T 

WATE2IT OOOl.AOO RELEASE FT05F001 ; DDEF FT05FOO 1 , VS , DUM . T5 , RET = T 

WATE2IT 0001500 RELEASE FT06F001 ; DDEF FT06FOC 1 , VS , DUM . T6 , RET=T 

WATE2IT 0001600 RELEASE FT05F001 ; DDEF FT08F00 1 , VS , DUM . T8 , RET=T 

WATE2IT 0001700 RELEASE FT09F001 ; DDEF FTO 9F0 0 1 , , BIN 

UATE2IT 00038C0 RELEASE FT13F001 ; DDEF FTl 3F0 0 1 , VS , S AVEFI L E 

WATE2IT 0001900 RELEASE FT1AF001;DDEF FTI AFOO 1 , VS , AMACPTS 

WATE2IT 0002000 IF ’ BOUT NONE ’ ; - 

WATE2IT 0002100 '"ElEASE FTl 0 FOO 1 ; DDEF FTIO FOO 1 , VS , ANSWER . BOUT , RET = $RET 

WATE2IT 0002200 IF ’ BOUT ’=* NONE’ ; RELEASE FTl OFOO 1 ; DDEF FTl OFO 01 , VS , ANSWER . BIN , RET=$RET 

WATE2IT 0002300 RELEASE FTl 5F00 1 : DDEF FT15F00 1 , VS , L ASTOPT 

WATE2IT 0002AOO RELEASE FTUF001;DDEF FTl 1 FOO 1 . VS , DUM . T1 1 , RET = T 

WAT12IT 0002500 RELEASE FT12F001 ; DDEF FTl 2F00 1 , VS , BMAPS 

WATE2IT 0002600 IF ’ BOUT *=* CONSOLE ' J REL EASE FTIOFOGI 

WATE2IT OOC2700 IF ’$ I H *-’ KEYBOARD R EL EASE FT09F001 

WATl'^IT 0002800 DDEF XXXI , VP , GRAPHICS , OPTION^JOBL IB 

WA7E2IT 0002850 IF ’ SB ACK * = ’ YES * : REL EASE GRAPHICS ; DDEF GRAPHICS, VS, DSNA«E=ZETA12C + 1 ) #DCB = (RECFM=V, LRECL=102B,MACRF=P) , RE* 
WATE2IT 0003000 IF LOADED= ’ TRUE GOTO ENGWT 
WATE2I7 00031 00 RELEASE X.YYY22 

WATE2IT 0003200 RELEASE DDNNEP;DDEF DDNNEP , VP , VCEl IB2 , OPTION= JOBLIB 
WATE2IT OOD3300 RELEASE WTNNEP;DDEF UTNNEP . VP , LIB . WATE2 , OPTION= JOBL IB 

WATE2IT 0003A00 IF * BHEWL I B * -= • NONE ’ ; REL EASE AANHEP;DDEF AANNEP,VP,$NEWLIB,OPTION=JOBLIB;D 'BNEWLIB LOADED* 

WATE2IT 0003500 LOAD BLOCKWT;D 'BLOCK DATA BIOCKWT IS LOADED' 

WATE2IT 0003600 QUALIFY ENGWT;D * ENGWT LOADED’ 

WATE2n 0003700 SET L 0 ADED= * TRUE’ 

WATE2IT 0003800 REMOVE AIL 
WATE2IT 0003900 QUALIFY NEPCAL 

WATE2IT OOOAOOO AT 1A0;IF NVOPT=0 & NIT>0;IF NIT=( NIT/5 )K5 ; D * ITERATING', NIT 
WATE2IT OOOAIOO AT 370C2);IF NVOPT=0 I NIT=51;D ’ B $ $ $ S FAILED TO CONVERGE $ $ B B B*,NIT 
WATE2IT 0009200 AT 570(2);IF NVOPT=0 I NIT<51;D •>>>>> SUCCESSFUL CONVERGENCE < < < < < *,HIT 
WATE2IT 0009300 AT 0(2);D ’ EXECUTION HAS BEGUN’ 

WATE2IT 0009900 DEBUG 
WATE2IT 0009500 QUALIFY WTE5T 

’WATE21T 0009600 AT 530;IF NVOPT=0;D WATENG,WATACC, MATE ' 

WATE2IT 0009625 QUALIFY BWTEST 

WAIEZIT 0009650 AT 30(2);Ir ’ BBACK ’ = * YES ’ ; SET SKIPIT= ' . TRUE . * 

WATE2IT 0CD9700 QUALIFY NEPCAL 
WATE2IT 0009800 ENGWT 
WATEZn 000990C CLOSE ANSWER. BOUT 
WATL2IT 00-350G0 CLOSE ANSWER.BIN 
WATE2IT 0005100 CLOSE LA5TOPT 
.'WATE2IT 00C5200 CLOSE BIN 
WATE2IT 0005300 CLOSE BMAP5 

WATE2IT 0005900 RELEASE FT;ERASE DUM.TljERASE DUN.T2;ERASE DUM.T8;ERASE DUM. Til 
WATE2IT 0005500 IF ’ BPRINT » = » YES ’ & ' BOUT ’-= ’NONE ’ ; PRINT ANSWER . BOUT ,,, EDIT 
WATE2IT 0005600 IF ' BPRINT - ’ YES ’ i ’ BOUT ’= ’NONE* ; PRINT ANSWER. BIN, ,, EDIT 
WATE237 0005700 IF ’ BBACK ’ = ’ YES ’ ; GTERM 
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puts to WATE-2 lor the engine to be redrawn and rewel^ed. These improvements would allow 
better use of the WATE-2 code in conjunction with NNEP. 
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